Negative specific heat in a quasi-2D generalized vorticity model 



T. D. AnderseCl and C. C. Lint 

Mathematical Sciences, RPI, 110 8th St., Troy, NY, 12180 
(Dated: February 1, 2008) 

Abstract 

Negative specific heat is a dramatic phenomenon where processes decrease in temperature when 
adding energy. It has been observed in gravo-thermal collapse of globular clusters. We now report 
finding this phenomenon in bundles of nearly parallel, periodic, single-sign generalized vortex 
filaments in the electron magnetohydrodynamic (EMH) model for the unbounded plane under 
strong magnetic confinement. We derive the specific heat using a steepest descent method and a 
mean field property. Our derivations show that as temperature increases, the overall size of the 
system increases exponentially and the energy drops. The implication of negative specific heat 
is a runaway reaction, resulting in a collapsing inner core surrounded by an expanding halo of 
filaments. 

PACS numbers: 47.27.jV, 47.32.cb, 52.25.Xz, 52.35.Ra 
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Negative specific heat is an unusual phenomenon first discovered in 1968 in microcanonical 
isolated system) statistical equilibrium models of gravo-thermal collapse in globular clusters 



13| . In gravo-thermal collapse, a disordered system of stars in isolation under-goes a process 



of core collapse with the following steps: (1) faster stars are lost to an outer halo where they 
slow down, (2) the loss of potential (gravitational) energy causes the core of stars to collapse 
inward some small amount, (3) the resulting collapse causes the stars in the core to speed up. 
If one considered the "temperature" of the cluster to be the average speed of the stars, this 
process has negative specific heat because a loss of energy results in an increase in overall 
temperature. Although one does not need to consider the cluster as a collection of molecules 
with negative specific heat to study its behavior, it can lead to insight into similar system, 
which has led us to this letter's findings. 

In the intervening four decades, negative specific heat has been observed in few other 
places. Its counterintuitive implication is that as a system loses energy its temperature 
increases and as one adds energy the system cools. In globular clusters, this results in 
a runaway reaction leading to a collapse of the cluster. In a magnetic fusion system or 
other thermally isolated plasma, should negative specific heat exist, the related runaway 
collapse could have profound implications for fusion where extreme confinement is critical 
to a sustained reaction. 

While 16| has proven that systems that are not isolated from the environment must 



have positive specific heat, the specific heat in isolated systems can be negative 12|. We 
now report finding negative specific heat in the thermally microcanonical system of confined 
generalized vorticity filaments, a finding of great interest to magnetohydro dynamics (MHD) 
and plasma physics. 

Our results have general applicability to vortex systems, and, although the calculations 
here-in presented have not yet appeared elsewhere, we have chosen a specific application to 
illustrate the value of our approach. In this letter, we address a plasma model known as the 
electron magnetohydrodynamical (EMH) model. 

Typically, magnetohydrodynamic plasma models are two fluid models, requiring equa- 
tions governing the electron motion and equations governing the ion motion coupled together 
[isl . The EMH model bypasses the two-fluid model by representing the electron fluid and 
the magnetic field as a single, generalized fluid with a neutralizing ion background that is 
stationary on the timescale chosen. 
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The EMH model takes the magnetic field, B = V x A, and the charged fiuid vorticity, 
u; = V X V, and combines them into a general vorticity field i7 = V x p where the generalized 
momentum, p = mv — eA, m is the electron mass, — e is the electron charge, v is the fiuid 
velocity field, and A is the magnetic vector potential field. For a brief overview of the model, 
see [l8|. A detailed model discussion can be found in [3i]. 

Our goal is to find the specific heat of this vortex model in statistical equilibrium given 
an appropriate definition for energy and a microcanonical probability distribution. Our 
approach is to describe the statistical behavior of a system of discrete, interacting vortex 
structures, rather than continuous vorticity, and consider the limiting case of an infinite 
number of infinitesimally weak vortex structures recovering the continuous vorticity field. 

The first step in our approach is to modify the continuous vorticity field, Q, such that 
it describes a large number, A^, of discrete vortex filaments, continuous in one dimension 
but point-like in the other two. Therefore, we assume that, instead of being continuously 
varying throughout, Q has a large number of periodic filaments that are nearly parallel to 
the z-axis. The nearly parallel assumption comes from the assumption of strong magnetic 
confinement and conservation of angular momentum which polarizes the filaments. The 
period, L, is assumed to be unity, L = 1, without loss of generality since all other distances 
and distance-dependent quantities can be scaled by L. 

The vorticity field, which depends on space, r e 3?^, now looks like: 

TV 1 

n{r)=J2 / rfrr,5(r-r,(r)), (1) 

^=l 

where r = (x, y, z) and rj = (xj, Zi). Periodicity requires that rj(0) = rj(l). Therefore, we 
assume that the discretization of r into has this property in Equation [H Because of the 
nearly parallel constraint, the arclength, r, has the property that r ~ Zi^li. For simplicity, 
we will represent rj(r) = (xj, r) as a complex number ipiir) = Xj(r) + iyiir). 

To simulate the vortex system on a computer, one must use a discrete field, but, even in 
our analytical approach, it is easier to start with a finite number of filaments and take the 
limit, N — i> oo, later, keeping total vorticity, A = Jsj^g f2(r)(ir, constant by rescaling. This 
is known as a non-extensive thermodynamic limit approach because the overall vortex 
strength stays constant even as the number of vortex filaments increase towards infinity. 
This common procedure is efficient in recovering the continuous vorticity field statistics. 

We use an approximate model for vortex behavior known as the local-induction ap- 



proximation (LIA), useful for nearly parallel filaments, combined with a two-dimensional 
logarithmic iiiteraction. Heuristically, this model comes from two separate results: One re- 
sult, that of [7], shows that a discrete field of perfectly parallel vortices has a logarithmic 
interaction. Another result for a single filament in 3D, that of [18], is a well-known first 
order approximation for the motion of fluid filaments, there-in extended to the EMH model. 
This LIA approximation, derived from the Biot-Savart law of magnetic/velocity induction 
for filaments of magnetism/ vorticity, depends on the binormal of the curve of the filament 
in 3-space causing Brownian variations. 



The combination of the two results of [3l and 



1^ yields the familiar London free energy 



of type-II superconductors [l5|, equally valid for generalized vorticity 18]: 

„1 N 

EN=a dry2T,/2\dUT)/dr\^- 

1 N N 

/ c?r5^^r,r,log|^,(r)-^,(r)|, (2) 

where a is the vortex elasticity in units of energy/length. The generalized angular momen- 
tum (angular momentum plus magnetic moment), from now on simply refered to as the 
angular momentum, is 

N 



J2^J dT\^p,{r)\'. (3) 
i=i -^0 



Fj is the vortex circulation. From this point on, we assume the vortex circulations are all 
scaled to unity, Fj = IVz, not entirely general but sufficient for the following results. 

Although we say the London free energy, we note that it is not describing the same type 
of system that it would describe in the type-II superconductor setting because these vortex 
filaments are generalized and not purely flux lines. 

What makes this system perhaps even more interesting is that it exists in an unbounded 
plane rather than with periodic boundary conditions, appropriate to superconductor research 
but not here, doing away with artificial planar periods. Other sources for this vorticity 
description as well as others can be found in 1J|, Q], [3, and 0, 11|. As far as we are 
aware, we are the first to apply it to the EMH model. 

For our isolated, classical system, the energy and angular momentum plus magnetic 
moment are conserved giving rise to the following probability distribution for the filaments 
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in equilibrium: 



P(s) = Z-^5(NHo - En -pMn)5(NR^ - Mn), 



(4) 



where Hq is the total "enthalpy" per vortex per period of the plasma, s is the complete state 
of the system, and Z — J ds5{NHo — En —pMn)5{NE? — Mn) is a normahzing factor called 
the partition function. Here En is the energy functional and Mn is the angular momentum. 
It is our intent to allow i?^, the size of the system, 

.1 



lim ( N-^ 



(5) 



to be determined by other parameters in the system and keep enthalpy and pressure, p, 
fixed. This distribution fully describes the density distribution of the filaments in the plane 
for every parameter. 

The size of the configuration space (partition function), Z, cannot be found in closed 
form by any known analytical methods. Since our aim is an explicit expression for specific 
heat, we need to know Z in closed form. 

To simplify the equations, we combine the large number of vortices into two average or 
"mean" vortices, which results in a mean vorticity field. This is the "mean-field" approach 
common in statistical mechanics. We also exploit the clear axisymmetry of the system and 
reduce the space from three to two dimensions. Given this approach, our mean vortices are 
as follows: one mean vortex is a mean distance from the origin. The other is the statistical 
center of charge of all the filamcnts-a single, perfectly straight filament fixed at the origin 
with strength of the remaining vortices, N — 1 ^ N. 

Given a filament i and a filament j, the mean-field approximation implies the following: 




where i is any filament index and the double bars, 
[0,1]. 

The energy function now reads as follows: 



dT\A{r)\' ^ m\, (6) 
1 1 , indicate £2-norm on the interval 



nl N 



a 




2 


dr 



(7) 



This assumption changes the interaction between vortices from fully coupled to fully de- 
coupled, that is, all vortices are statistically independent, and the statistics of the entire 



system can be found from the statistical behavior of a single vortex. We modify H] and [3] 
appropriately and, from here on forth, drop primes since neither the full energy nor the 
complete distribution nor the full angular momentum will be further referenced. 

In statistical mechanics of isolated systems, all equilibrium statistics can be determined 
from maximizing the entropy. The entropy per filament, Sn, is defined by. 



j Dij6{NHo -En- pNR^)5{NR^ - Mn). (8) 



This definition implies, 

Sn = log 



j DipS{NHo -En- pNR^)5{NR^ - Mn) 



(9) 



We have now set the stage to describe our derivation of negative specific heat. 

Given the space available, we proceed to outline, rather than fully derive, our method 
of obtaining an explicit formula for the maximal entropy of this mean-field system in the 
non-extensive thermodynamic limit (as defined above) from which we obtain an explicit, 
closed form formula for the specific heat. A full description of the mathematics involved is 
forthcoming in another paper. 



First it is important to note that our approach re 



ies heavily on the steepest descent 



llTj . The works of 9|, |lO|] have pre 



methods in p] and spherical model approach in 
ceded and inspired this work in their novel applications of the spherical model to barotropic 
vorticity models on the sphere. These works lay the ground for our derivation. 

The steepest-descent and consequently spherical model methods convert Dirac-delta func- 
tions into their Fourier-space equivalent integral representations. The advantage of this inte- 
gral representation is that one can perform steepest-descent evaluation on it. Therefore, the 
procedure for our derivation relies on two closely related facts: the integral representation 
of the Dirac-delta function in a microcanonical distribution, for example: 

6{Nx - Nxo) = / |^eAf/3(-xo)^ (^q) 
J/3o-ioo 2vr2 

and the steepest-descent limit, again only an example: 

lim — / d/?e^^(^-^°) = e*(^-^'»), (11) 

which allows quantities such as entropy to be determined, provided we can determine what 
(3q is, a method for determining which provides. 



Using the first fact, we can convert tlie microcanonical problem into tlie canonical prob- 
lem by converting delta functions into integrals and re-ordering the phase-space {tp) and 
parameter-space (/?) integrals to arrive at 



e^- = / ^e^'^^^Z,^,,, (12) 
where 

Zbath = J Dtlje-^^^'''''^5iNR^ - Mn) (13) 

is a canonical-like partition function for the system in an external heat-bath but with a 
microcanonical angular momentum constraint. In the infinite N limit, Zbath approaches the 
canonical partition function (i.e. the same function but without the Dirac-delta factor in 
the integrand of Equation [T^ . making their use interchangeable in the limit. 

To do this, i.e. convert the microcanonical problem to a canonical one, one has to prove 
that the integrals are finite, and indeed we can but do not show it here. One this is done, 
solving for the canonical partition function becomes the solution for the microcanonical 
partition function in steepest-descent, and, as is well-known, exponentials are much easier 
to integrate than Delta-functions over complicated energy surfaces. 

The maximal entropy per filament per period, having the form, 

SmaAHo) = lim N-'Sn, (14) 

is where the limit comes into play. 

One problem commonly seen with logarithmic. Coulomb-style interactions such as this is 
that as the number of "particles" (in this case vortex filaments) grows infinitely large, the 
interaction energy grows with the square of the number of filaments. This type of super- 
extension is never seen in traditional molecular system where there is a clear distance cut-off, 
but it is seen in both Coulomb and gravitational systems like vorticity models and globular 
clusters. The solution is to rescale the temperature, which in turn, rescales the interaction 
energy. Rescaling the temperature causes a chain of necessary scalings to restore the balance 
so that other quantities do not go to zero: /?' = PN, a' = a/N, p' = p/N, and Hq = Hq/N. 
These are, of course, reasonable because only the interaction energy needs a rescaling. 

With all these scalings there are no more mathematical obstructions, and the maximal 
entropy can be found by the standard mathematical procedures in Q and jl|. We provide 



only the final formula obtained in view of space constraints, but the procedure is quite 
straightforward once the appropriate framework is set up: 

SmaAHo) = [^M + f log(i?') - - f^oP'R'- (15) 

where 

^ = 8^' ' (^^^ 

where the mean temperature, Pq, is as yet unknown. This entropy is exact within the 
mean-field assumption for N —>■ oo. 

By jsl, we find the unknown multiplier, /3o, by relating the enthalpy per filament param- 
eter, Hq, to the mean enthalpy, NHq = {E^ +pMn), where (■) denotes average against 
Equation HI 

By Equation m the average enthalpy is given by 

+ ^^^^ - / D^6iNHo -Ej,- pMj,)6iNR^ - M^) ' ^^^^ 



and, again going through some steepest-descent-based calculations given in we find the 
formula, 

exactly. We cannot give an explicit expression for /5o because it is a root of a transcendental 
equation, but such is unnecessary for the following negative specific heat result: 
We define specific heat at constant generalized pressure, p, 

and after evaluating with Equation [18] and simplifying (dropping primes and 0-subscripts) 

c, = gf, -iV (20) 

Equation [20] is significant. It indicates that the specific heat is not only negative for this 
system, but strictly negative if parameters are non-zero (Figure [1]). In the low-temperature 
(large f3) case, for constant field strength, does not change significantly with temperature 
indicating that filaments are in a stable configuration for a large range of low-temperatures. 
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FIG. 1: The specific heat at constant pressure (Equation I20p for the thermally isolated system is 
negative, meaning that the constant pressure Enthalpy Per Length (Equation [T8|) decreases with 
increasing temperature. (Here a' = 5 x 10^ and p' = 8 x 10^.) 

Because the filaments do not move relative to one another at low-temperatures and the 
self-induction is negligible, the enthalpy does not change. When temperature becomes high 
the internal entropy causes a massive expansion in the overall size of the system (Figure [2]), 
and energy of the logarithmic interaction decreases far more than the enthalpy of the self- 
induction increases. The strong magnetic field absorbs this energy, but, since it is assumed 
to be an infinitely massive reservoir able to maintain the enthalpy at Hq, the confinement 
remains constant. 

On a side note, because the specific heat in Equation [20] does not cross the axis (i.e. is 
never positive) for any positive parameters, the expansion in in the canonical system is 
not likely a phase transition in /?, but a continuous "transition" jl2|. The system appears 
to change behavior significantly, but it changes without any discontinuity in the free energy 
or specific heat. 

That only a local-induction approximation is necessary for this analysis is interesting 
because it is so simple and so moderately different from 2D vortex system which do not 
exhibit negative specific heat in statistical equilibrium. What happens in the nearly parallel 
vortex filament case can be explained in steps similar to the catastrophic gravo-thermal 
collapse: (1) a vortex's Brownian motion causes it or part of it to move away from the 
center, (2) potential energy decreases, (3) the vortices in the center can move closer together 
and temperature increases. 

As mentioned in the preceding paragraph, the negative specific heat indicates a runaway 
reaction (i.e. the fixed energy, fixed angular momentum equilibrium point is meta-stable) 
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FIG. 2: The Mean Square Vortex Position (Equation I16p increases exponentially at high temper- 
ature, while it is nearly constant at low-temperature. (Here a = 5 x 10^ and p = 8 x 10^.) 

in the microcanonical system that models strongly confined, single-sign generalized vortex 
filaments. This is the insight we have from gravo-thermal collapse: We hypothesize that 
this kind of meta-stability, observed in gravo-thermal collapse of globular clusters [l^, can 
lead to two possible outcomes: (1) a collapse similar to globular clusters in which an outer 
halo of columns separates from an inner core that collapses in on itself, possibly resulting in 
nuclear fusion, (2) a turbulent expansion of the entire system. Further research will focus 
on answering this question, but clearly 3D effects, even only a LIA, are crucial. 
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